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Abstract 

We establish necessary and sufhcient conditions for convergence (in the sense of finite dimensional 
distributions) of multiplicative measures on the set of partitions. We show that this convergence is 
equivalent to asymptotic independence of finite sizes of components. The multiplicative measures 
depict component spectra of random structures, the equilibrium of classic models of statistical me- 
chanics and stochastic processes of coagulation-fragmentation. We then apply Schur's tauberian 
lemma and some results from additive number theory and enumerative combinatorics, in order to 
verify the conditions derived in important special cases. Our results demonstrate that the common 
belief that interacting groups in mean field models become independent as the number of particles 
goes to infinity, is not true in general. 
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1 Probabilistic setting and its applications 



Let {Zj, j > 1} be a sequence of independent integer valued random variables that induces a sequence 
of random vectors {K*-"^ = {K^\ . . . , K^'^), n > 1} given by 

n 

£(k(")) =£(Zi,...,Z„|^jZ,- = n), n = l,2,.... (1.1) 
It follows from ([LT]) that K^") G r2„, n > 1, where 

n 

^n = {1] = ih, ...,kn) : ^ jfcj = n} (1.2) 

is the set of all partitions rj of an integer n. In probabilistic combinatorics, (jl.ip is called the condi- 
tioning relation (see[3]), while the sequence of vectors {K^"), n > 1}, is called the counting process. 
Next, denote by the probability measure on Qn induced by the conditioning relation (jl.ip : 

finiv) ■■= ]P(K(") = r/), r]£ n„, n > 1 

and let 

ajf'^ = TP{Zj = k), k> 0, j > 1. (1.3) 



We then have 



where 



n 



/i„(7/) = r? = (/ci, . . . ,/c„) G f]„, n > 1, 



cn = FiY^jZj = n) = E n % ' V = iki,...,kn) (1.4) 

J = l »?Gn„ jr = l 

is the partition function for the measure /U„. We will assume throughout this paper that the proba- 

(7) 

bilities , k >0, j > 1 are such that c„ > 0, n > 1. 

It is clear that the sequence of measures {fin, > 1} induced by (jl.ip is uniquely defined by the 
array of probabilities {0'^^\ J ^ 1; k > 0}. However, this correspondence is not a bijection. In fact, 
the "tilting" transformation (see [3]) of the probabilities: 

"' ''" = Safe' ' ^ ^ 

where p > and S^^\p) is the normalizing constant, does not change the sequence {pn n ^ !}• But 
this transformation does affect the partition function c„ : 
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Note that the tilting is defined for all finite p > 0, such that 

oo 

50)(p) = ^^'=4^)<oo, j>l. (1.7) 

fc=0 

It is a remarkable fact that the representation (jl.ip provides a mathematical formalism for a vari- 
ety of models in seemingly unrelated contexts. Let us briefly describe four main fields of application 
of this setting. 

• Decomposable combinatorial structures (for more details see [3] , [1] , [27] ) . We consider 
combinatorial structures. The size of such a structure is defined to be the number of elements in 
it. A decomposable structure of size n is a union of indecomposable components, so that the counts 
ki, . . . ,kn of components of sizes 1, . . . , n respectively, form an integer partition of n. It is assumed 
that each component of size j belongs to one of rrij types. The three classes of decomposable struc- 
tures: assemblies, multisets and selections, encompass the whole universe of classical combinatorial 
objects. Assemblies are structures composed of labelled elements. The class of assemblies includes 
permutations decomposed into cycles {rrij = (j — 1)!), forests composed of rooted trees with labelled 
vertices {nij = ), graphs composed of connected subgraphs with labelled vertices (m, ~ 2(2)), 
etc. Multisets are formed from unlabelled elements. Examples of this class of structures are integer 
partitions (mj = 1), planar partitions {nij = j) and mapping patterns {nij ~ p = 0.3383). Fi- 
nally, selections are defined as multisets with distinct components, which means that all component 
counts fcj, j = 1, . . . , n are either or 1. A typical example of a selection is an integer partition into 
distinct parts {rrij = 1). 

The basic problem in enumerative combinatorics is to find the asymptotics, as n goes to infinity, of 
the number of a certain class of structures of size n, having a given component spectrum {ki , . . . ,kn) £ 

As a part of this question, the asymptotics of the total number of given structures of size n is of 
special interest. 

The starting point of the probabilistic method considered is the definition of a random structure 
of size n which is a random element n„ distributed uniformly on the finite set of all given structures 
with size n. Next it is defined the induced random component spectrum (also called counting process) 

K(-) = K(")(nO = (Ki"\ . . . , K(")), n > 1, 

(n) 

where the random variable represents the number of components of size j in n„. 

It turns out that the representation p.l|) of the distribution of K^") is valid for the aforementioned 
three classes of combinatorial structures. Namely, assemblies, multisets and selections are induced 
respectively, by the following three types of random variables Zj, j > 1 : Poisson {Po{aj), aj = -jf ), 
Negative Binomial {N B{pj,mj)) and Binomial {Bi{pj,mj)). 
• Models of ideal gas (for references see [33], [25], Ch.l2, [31]). 
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An ideal gas, treated as a quantum system, is a collection of perfectly elastic particles (atoms or 
molecules) which collide but otherwise do not interact with each other. It is assumed that the total 
internal energy E of a gas is the microscopic energy of the random motion of individual particles and 
that E is partitioned between the particles, so that kj, called an occupation number, is the number 
of particles with energy value (energy eigenvalue) ej, i.e. 

E = Y,^jkj. (1.8) 

In the case of a d- dimensional ideal gas, a particle is viewed as a lattice point and the energy 

values are of the following special form: 

d 

e||q||2 = c||qf, q = (/i, . . . , y e Z^, ||qf = ^ ll (1.9) 

s=l 

where c > is a known constant that does not depend on q. 

Now we see from ()1.8p . ()1.9p that, given a total energy E, a state of an ideal gas is determined by 
a partition t] = (ki, . . . , kn) of an integer n. The probability distribution of the energy states rj which 
varies from model to model, is defined by a measure on the state space Qn- By laws of statistical 
mechanics, these measures are forced to be of the form (jl.4p . with the numbers defining the 
type of a model of ideal gas considered. Vershik ([33j) suggested to call the class of measures ()1.4p 
multiplicative. Observe that the multiplicative form (jl.4p of the sequence of measures fin is implied 
by the fact that the random variables Zj, j > 1 in (|l.ip do not depend on n. 

The following three basic models (statistics) of an ideal gas are distinguished: 
Maxwell-Boltzmann(MB) (labelled particles, without restrictions on values of occupation numbers), 
Bose-Einstein(BE) (indistinguishable particles, without restrictions on values of occupation numbers) 
and Fermi-Dirac(FD) (indistinguishable particles, such that no more than one particle may occupy 
a given energy level). In full similarity to the setting for combinatorial structures, MB, BE and FD 
statistics are induced by the same three types of random variables Zj in (jl.ip : Po{aj), NB{pj,nij) 
and Bi{pj,mj) correspondingly. In the case of the classic 2- dimensional MB, BE and FD statistics, 
the parameters nij = const, j > 1, while to the parameters pj and aj it is attributed the meaning 
of the average of the occupation number at energy level j which has a known in quantum mechanics 
expression via the corresponding energy eigenvalue ej. This says that the classic BE and FD models 
conform respectively to geometric and bernoulli distributions of Zj, j > 1, which induce the uniform 
probability measures fin on the set of all integer partitions of n and on the set of integer partitions 
of n with distinct summands, respectively. 

In the case of a 2 < dimensional ideal gas, the situation is more complicated. By (jl.Sp . (jl.Op . 
to each energy level j is naturally prescribed a "weight" r^ij) which is the number of representations 
of the natural number j as a sum of d integer squares. In other words, rd{j) is the number of distinct 
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lattice positions q = {li, . . . ,1^) S Z"^ of a particle on a sphere of radius ^/J, i.e. with the energy 
level j. It is known from number theory (see e.g. [24J) that for d > 5, 

Cij"^ <rd{j) < C2j^, j > 1, 

where Ci,C2 are positive constants depending on d, and that for d = 2,3,4 the functions rdij) 
oscillate wildly (in j), while, obviously, 

{2, if 7 is a square 
(1.10) 
0, otherwise. 

Employing known properties of r^lj), an important fact was proven in [31] that for the sake of 
asymptotic analysis, it is possible to treat the d- dimensional BE and FD models as the classic ones 
with parameters rrij = cj^ , (3 > —1, c > 0, where /? = if d > 2 and /? = — ^, if d = 1. 

• Coagulation fragmentation processes on the set of integer patitions (see [I5].[3]). 
We will show that a multiplicative measure can be viewed as an equilibrium of a classic coagulation- 
fragmentation process (CFP) which is a time-continuous Markov chain on defined as follows. A 
state ry = (ki, . . . , kn) G of a CFP depicts a partition of a total number n of identical particles 
(animals, atoms, stars, human beings, etc) into clusters (groups) of different sizes, so that kj is the 
number of clusters of size j. The only possible infinitesimal (in time) transitions are coagulation 
(merging) of two clusters of sizes i and j into one cluster of size i+ j and fragmentation (splitting) of 
a cluster of size i+ j into two clusters of sizes i and j. Given a state rj S r2„ with ki ,kj>0 for some 
^ ^ i, j ^ n, denote by r/^*'-'^ G 0„ the state that is obtained from r] by the coagulation of any two 
clusters of sizes i and j, and denote by Ucirj,"!]^^'^^) the rate of the infinitesimal transition ij rj^^'i) , 
Similarly, for a given state rj G J7„ with ki^j > 0, let 'r]{i,j) be the state that is obtained from r/ by 
the fragmentation of any cluster of size i + j into two clusters of sizes i and j, and let Uf{rj, Viij)) be 
the rate of the infinitesimal transition rj — > ^/(jj). Denoting by 

the ratio of the above transitions, the important property of reversibility of multiplicative measures 
is derived by verifying the detailed balance condition. 



Proposition 1 A multiplicative measure fin defined by is reversible with respect to the transi- 
tion rates Uc,Uf, such that their ratio satisfies: 



(i) (2i) 
1 n 
X-2%» + l 
(i) (2i) : 



if j ■ ki-,kj > 



(1.11) 
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An immediate consequence of the Proposition 1 is that a multiphcative measure fin defined by ()1.4p 
is the equihbrium distribution of a CFP with transition rates obeying the condition (jl.lip . 
We now distinguish a class of CFP's with transition rates of the form 

Uf{v,V{i,j)) 

where (j), (/? are some symmetric nonnegative functions on the set of pairs of positive integers. Treating 
the functions (j), ip in (|1.12p as the rates of a single coagulation and a single fragmentation respectively, 
the induced CFP's can be viewed as mean-field models on the set In fact, (jl.l2p tells us that 
at any state rj £ Qn, each cluster can coagulate with each other one or can be fragmented into two 
parts, so that the net rates of the transitions r/ — > ^(j j) and rj 7y(*j) are sums of the rates of all 
possible, at the state rj, single coagulations and fragmentations respectively. It is proven in [T2] with 
the help of Kolmogorov cycle condition that in the case when all rates of single coagulations and 
fragmentations are positive, the CFP's given by (|1.12p are reversible if and only if the ratios of the 
single transitions are of the form 

^ = ^, ^,J>1, (1.13) 

with some Oj > 0, i > 1. The corresponding CFP's are known as classical reversible models of 
clustering and networks studied in 1970-s by Kelly and Whittle (see |15j-|17j and references therein). 
The equilibrium measures of the mean field CFP's with the rates ()1.12p . (jl.lSp are multiplicative 
measures /x„ induced by the conditioning relation (jl.ip with Zj distributed Po{aj), j > 1. A clear 
example of a reversible CFP which is not a mean field, is provided by letting Zj = Gj — 1, where Gj 
is distributed geometrically {p^ , < p < 1). Then a^j^^ = p^^qj, Qj = ^ — p' , J > 1, k > 0, and the 
corresponding measure /i^ is a uniform one on the set while reversibility holds iff q{i];i,j) = 1, 
by (jl.lip . Simple calculations demonstrate that, due to the last fact, the representation (|1.12p of 
net transition rates fails for the CFP considered. 

• CFP's on set part it ions ([7], [50]. [S]). We assume here that in the preceding set up for CFP's, 
particles are labelled by 1, . . . ,n, so that the state space of the system of clusters related to some 
CFP, becomes the set = {vr[„]} of all partitions -7r[„] of the set [n] = {1, . . . ,n}. Recall that a 
partition of [n] into k blocks (clusters) Ai, . . . , is iT[n],k = (^i; • • • > ^k), where Aj, 1 < j < k < n 
are nonempty and disjoint subsets of [n] whose union is [n] and which are numbered, e.g. in the 
order of their least element. Denoting \Aj\ the size of a cluster Aj, we further assign to each Aj, a 
weight m^Ajl which is a number of possible states of Aj, the states can be e.g., shapes (in the plane or 
in space), colors, energy levels, etc. This says that to the set partition 7r[„] ^ correspond 0^=1 "^|Aj 
different structures with the same blocks Ai, . . . , Ai^, so that the total number of structures formed 



kikj<f>{i,j), if i^j, kikj > 

ki{ki - l)4>{i,i), if i = j, ki > 2, 



ki+j(p{i,j), 2 < i + j < n, ki+j > 2, 



(1.12) 
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by all partitions of the set [n] into k clusters is equal to 

k 

Yl llrn\A^\ ■■= Bn,k, (1.14) 

where Bn^k is known as a Bell polynomial in weights mi, . . . ,mn-k+i- Similar to the above setting 
for decomposable combinatorial structures, a random structure Hj^]^^ is the one chosen randomly 
from the set of Bn,k structures. As a result, for a given k a measure P[n],k on the set ^[n],k = {''^[n],k} 
is induced: 

. . ni=i"^|A,| 

P[n],k['^[n]) = ^ , 7r[„] G (1.15) 

In a more general setting which encompasses a variety of models (see [7].|30)). the weights rrij in (jl.lSp 
are allowed to be arbitrary nonnegative numbers. Pitman [30] calls the Hf^] a Gibbs partition and 
the measure fc microcanonical Gibbs distribution. Obviously, the vector (|Ai|, . . . |Afc|) of block 
size counts defines a partition of the integer n into k summands induced by the generic set partition 
TT[n],kj aiid it is known that to each t] = {ki, . . . , kn) G f^n, such that /ci + . . . + /c„ = /c, correspond 

n! 

n-=i(%o(j!)'^ 

different partitions ir^^.k S ^[n],k^ each one of them having the same probability P[n],fc(^[n],fc)i given 
by (jl.lSp . Thus, the Gibbs distribution on ^[n]^k induces the Gibbs distribution p„ ^ on the set 
^n,k of integer partitions of n into k positive summands: 

" m 1 

p^^,,{ri) = {Bn,kr^X{C^f^jTrT-v rj = {ki,...,kn)enn,k, (1-16) 
where the partition function -Bn.fc defined as in (jl.l4p can be rewritten in the following form: 

n 

From (jl.l6p it is easy to derive Kolchin's representation of Gibbs partitions (see [30], Theorem 
1.2). From the other hand, the distribution pn^k given by (11.161) is produced by conditioning the 
multiplicative measure measure /i„ given by (jl.3p . on the event Zi + . . . + Zn = k, with Zj ~ 
Po{aj), aj = j > 1. However, this embedding of the generic model associated with set partitions 
of [n] into the setting for conditioning relation (|l.ip does not facilitate the study of a wealth of 
problems (see |7j) arising from treating P[n],fe, A: = 1, . . . ,n as marginal distributions of irreversible 
time continuous markov processes of pure fragmentation (or pure coagulation) on the state space 
[n]. The study of these problems was initiated by Kingman and Pitman and is being extensively 
continued by a group of researchers including Pitman, Bertoin, Berestycki, Gnedin et al. 

In what follows we will refer to all models induced by the conditioning relation (jl.ip as random 
structures. 
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2 Objective and Summary. 

In this paper, we study the asymptotic behaviour of the random vector {k[^\ . . . , K^^^) of smaU 
(to compare with n) component counts, defined by as n — > oo and Z > 1 is fixed. Based 

on the independence of the random variables Zj, j > 1 in (jl.ip . it was a common behef that the 
counts become independent, as n — > oo. We show that this is not true in general. Our main result 
which is Theorem 1 in Section 3, consists of proving that the asymptotic independence of small 
component counts is equivalent to the weak convergence of the above sequence of random vectors, 
and of establishing a characterization of convergent models (multiplicative measures). This is done 
by employing the theory of RTp sequences which originated from Compton's work on logical limit 
laws and has been further developed by Burris and Bell. Based on Theorem 1, Schur's tauberian 
lemma and results on asymptotics of random structures, we answer in Section 4 the question of 
convergence of counting processes for the three types of basic models distinguished in Section 1. 
It turns out that many models related to multiplicative measures are divergent. In a parallel way 
we discuss the problem of convergence for CFP's. The last Section 5 contains concluding remarks, 
among them a historical background of the discussed problem. 

3 Main result. 

Definition 1 . We say that the counting process {K*^"\ n > 1}, is convergent, if for each finite 
I > 1, the probability law C{k[^\ . . . ,K^"'^) weakly converges, as n ^ oo, to some probability law Fi 
on M', / > 1. Moreover, we say that counts k["\ . . . jK^"'^ of small components of the random vector 
Qj^g asymptotically independent if the above laws Fi are product measures on M', for all finite 

I > 2. 

Note that in contrast to the setting for limit shapes (see e.g. [33], [l8]), in this paper we study the 

weak convergence of non scaled multiplicative measures. 

For a fixed / > 1, given ki, . . . ,ki and sufficiently large n, we denote 

I n-Mi 

Ml = ^jkj, and T^'\n, Mi) = P( jZ, =n- Mi). (3.18) 

j=i j=i+i 



It is immediate that 



r«(n,MO =T«(n-M,,0) :=T^'l^,, / > 1, > 0. (3.19) 



^n-Mi ' 

Assuming in what follows that 

ai^ > 0, j > 1, 
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we will be dealing with the "scaled" quantities d^^\ Cn, and defined by 



a? = \, k>0, j>l, (3.20) 



c„ = ( n ao^) n > 1, CO = 1, (3.21) 



3=1 

n—k 



j=i+i 

It will become clear from our discussion in Section 4 that in the case of classic decomposable 
structures, the quantities c„ and T^\^ have a significant combinatorial meaning. Namely, p~^Cn is 
the number of structures of size n in the case of multisets and selections, whereas n\cn is the number 
of assemblies of size n. In analogous way, T^'l^ is related to the number of structures of size n — k 
with all component sizes greater than I. 

With the help of the above notations, we have 

P(<) = h,..., K'i"^ = k) = ( n <>( E JZ, =n-M,) = 

j=i j=i+i 

I fii) 

where in the last step we have used the fact that 

n Ml 

P( E jZj =n-Mi)= T^l^,^ n ^iZn-M,+k = 0) 

j=l+l k=l 

and the definitions (|3.20p . (j3.2ip and (j3.22p of the scaled quantities. Note that in view of (|3.18p . 
'^n-Mi '^^^^ al\ki,...,ki: Y!j=i3kj = M/. 

Central to our subsequent study is the notion of smoothly growing real sequences RTp which 
definition we adopt from [lOj . [6]. 

Definition 2 . i?Tp, Q < p < co is the collection of sequences {dn}n>i of nonnegative numbers 
that satisfy 

lim = p. (3.24) 

n-^oo dn+1 

m 

Sequences in RTp play a key role in Compton's theory of logical limit laws and in additive number 

theory (for references see [TO], [5], [6]). 

Now we are prepared to state our main result. 
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Theorem 1 . The conditions (a),(h) and (c) below are equivalent: 
(a) The counting process {K^"\ n > 1} is convergent; 

(h) Counts of small components of random vector K.^"^^ are asymptotically independent; 
(c) {cri}n>o £ RTp^ for somc < p < oo, and 



rfil) I - 1 



< g(') := lim ^ = ( TT S^^\p)] < oo, / > 1, (3.25) 

n^oo Cn \ J 



Where ~S(^){p)=Y.T=oP^''~4^- 



Proof. 

Assumption (a) implies that the fraction in the RHS of (|3.23|) has finite Hmits, as n ^ oo, for all 
fixed Ml > 0. Moreover, for any / > 1 there exists a Mi > 0, such that 

f(i) 

lim > 0. (3.26) 



We write now 

-^n-Mi-l _ I -^n-Mi-l \ ( C„_i '\ 



Cn-1 



f.(0 

n 



to conclude from the above two facts that Cn E RTp^ for some ^ p ^ oo. Consequently, lim^ — — ; — 
g^') < oo, / > 1, and we get 

f{i) 

lim ^^=g«/^', />!, (3.28) 

n— >oo Cn 

for all < Ml < oo. As a result, < g*-'^ < oo and the limit, as n ^ oo of the RHS of (j3.23p equals 
to ^ 

q^'^fld^ll^fP"^. (3.29) 

(j3.29p together with (j3.23p show that the probability limit law, as n ^ oo, of the vector of small 
counts (ifi"\...,Kj"^) is a product measure on W, for all / > 2. 

Thus, we proved the implications (a) =^ (6) and (a) =^ (c). Now the proof of the converse of the 
second of the two implications is obvious. ■ 

Remark 1 (i) In condition (c), the case p = is not excluded. Setting 0*^ = 1 and recalling that 
clq ^ = 1, j > 1, it follows from 113. 23]) and 113. 29\) that in this case the limit law of {k[^\ . . . , -fc/"^), 
for I > 1 is the measure concentrated at the singleton (0, . . . , 0) G i?', whereas q^^^ = 1, in accordance 



T. 



(0 



with (MW> and lim„^oo = 0, for all Mi > 



(a) Condition (c) also implies {Ti'^}n>o £ RTp, for all I > 1, with the same < p < oo as for the 
sequence {c„}„>o- This can be seen by writing 



ft ^ lfl_\ (it, 
~^ \ftJ \ , 
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and then applying the fact that < < oo and (E 



(in) One can see from the proof of Theorem 1 that condition g^'-* > 0, I > 1 which is a part of the 
condition (c), ensures the tightness of the corresponding sequences of finite dimensional probability 
measures. 

To formulate the forthcoming corollary, we need to extend the definition (jl.Sp of the tilting 
transformation to the case /? = 0, in the following natural way: 

a?(0)= (3.30) 
I 0, otherwise. 

Corollary 1 Let K^") be a convergent counting process, such that Cn € RTp, for some < p < oo. 
Then 



lim P(irJ") = ki,..., k/"^ = ki) = ll 4f (/5), / > 1, (3.31) 



kj 

Ji) 



where (p) are the generic probabilities tilted with the p. 

Proof. By ([3:25]) . (13:29]) and the definitions (fT3 ]) . ([3^0]) and dS^Oj), it follows from Theorem 1 that 
for a convergent counting process, 

hm P(<) = fci, . . . , K'f'^ = h) = T\ - 



^^S(j){p) 



I cP^ ' 



Remark 2 (i) Corollary 1 says that asymptotic independence with the limit product measure com- 

(i) 

posed of generic probabilities takes place only if p = 1; 

(a) We denote by Cn{0) the quantity Cn corresponding to the tilting of the probabilities a^^'^ with 
6 > and recall that a multiplicative measure pn is invariant under all possible tiltings of the 
probabilities with > 0. By definitions d j.5]) and 13.21\) we then have 

^) = ^"c„, n > 0. 

So, if {cn, n > 0} G RTp, for some < p < oo, then 

{cje), n > 0} G RTe. 

This clarifies the meaning of Corollary 1: Civen a convergent generic counting process such that 
{cn, n > 0} G RTp, with some p > 0, the whole family of counting processes obtained by tilting the 
original one with all possible 9 > 0, has the same limit finite dimensional distributions as the counting 
process tilted with the above p > 0, so that the corresponding quantity {cn{p), n > 0} G RTi. 
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4 Convergent and divergent random structures. 



We agree to call a random structure convergent (divergent) if the corresponding to it counting process 
converges (diverges) respectively, in the sense of Definition 1. 

Our tool for verification condition (j3.25p of Theorem 1 for the models considered, will be the remark- 
able Schur's tauberian lemma cited below. With an obvious abuse of notation, we say that a power 
series f{x) = J2n>o '^n^" is in RTp if {dn}n>o G RTp. We denote by * the Cauchy product, which is 
a usual product extended to formal power series (see |10j). 

Lemma 1 (Schur (1918), see lWI,p. 62). 

Let f = fi* f2, where f, fi, f2 are power series with coefficients dn, dn\dn\ ''^ > 0, respectively, 
such that : 

(a) fi £ RTp for some < p < oo and 

(b) the radius of convergence 0//2 is greater than p. 
Then 

hm 4) = (4-32) 

We also note the fact that /i G RTp together with (j4.32p imply / G RTp. This can be seen by 
writing 

dn—l f "n— 1 \ / "n— 1 \ f '-I'n 



Let S^^^ be the generating probability function of the random variable jZj in the conditioning relation 



fc>0 

Clearly, the radius of converge of S^^^ is > 1, for all j > 1. 

It follows from (jl.4p that the Cauchy products g := Ylj>i S^^^ and gj^(i) := nj>/+i ^^''^ 
generating functions for the sequences {cn}n>o and {Tn^}n>o, I ^ li respectively. This fact is a 
characteristic feature of multiplicative measures on the set of partitions. In accordance with the 
preceding notations, we also denote by S^^\ g and gf^i) the generating functions for the scaled 
sequences {al^\ ^ > 1}, {cn, n > 0} and {Tn\ n > 0} respectively. Finally, denoting g^''^ = 
ltj=iS^^\ / > 1 we have 

9 = m) * 9^'^ (4-33) 

or, equivalently, 

1 

5( 



~9fin = ~9*{^ ■ (4.34) 
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Supposing that condition (c) of Theorem 1 holds and employing (ii) of Remark 1, allows to apply 
Schur's lemma either to ()4.33p or to (j4.34p . depending on whether g^''^ or have the radius of 
convergence greater than p. As a result we find or q^''^ respectively for the random structures 
considered. 

In what follows we examine the convergence of counting processes for the three basic types of 
random structures: assemblies, multisets and selections, described in Section 1. Furthermore, the 
results obtained explain the crucial difference in the asymptotic behaviour of mean field and non 
mean field CFP's associated with the above structures. 

The implementation of Theorem 1 for random structures with given sequences of parameters 
{n-j) j ^ 1} requires validation of the condition g £ RTp for some < p < oo. Our treatment of 
the problem is based on known sufficient conditions on sequences {raj, j > 1} providing the above 
RTp property for the corresponding sequences {cn, n > 0}. The conditions we employ are products 
of two quite different lines of research. The first is Burris-Bell theory ([5], [6]) of RTp sequences, 
developed with the help of analytical tools stemming from Tauberian theory. Motivation of this 
research came from Compton's (1980s) theory of logical limit laws and also from the additive num- 
ber system theory. The second source of the sufficient conditions are asymptotic formulae for the 
number of decomposable structures, in particular recent results by Barbour, Freiman, Stark and the 
author, derived with the help of probabilistic methods. Each one of the two approaches has its own 
limitations, some of which are noted later on in this Section. In this connection, we observe that 
there is little hope to obtain plausible necessary conditions for the RTp property of g. 
The above mentioned sufficient conditions combined with Theorem 1, allow us to obtain character- 
izations of convergence for each type of random structures, and to provide examples of convergent 
and divergent structures. 

Firstly, following [3| we call regularly varying the random structures induced by random variables 
Zj, j >1 in ([LTD, with 

TEZj ~ const y^f, a G R, y > 0, j ^ oo. (4.35) 

Since the asymptotic behaviour of regularly varying random structures appears to be in accordance 
with the behaviour of the series X^jli (s^e [S], [TO], [1], [23]), it was suggested in |4] to distinguish 
the following three types of regularly varying structures: logarithmic (a = —1), convergent (a < — 1) 
and expansive (a > —1). (It goes without saying that in this classification, the meaning of a 
convergent structure is different from the one in the present paper). As in [TO], we extend the above 
definition of expansive structures to include random structures with EZj, j > 1 oscillating (in j) 
between two regularly varying functions, namely 

(EZ,, j > 1) G Hn,r2;y) := {/ = f{j) ■■ iiy'p-^ < f{j) < i2y'p-\ j > i, y > 0}, 
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where 7^, i = 1,2 are positive constants and < ri < r2. (The requirement ri,r2 > is the 
characteristic feature of the expansive case). 

Assuming the conditions of Schur's lemma are fulfilled, we then get 



^(0 



{g^Hp)y\ (4.36) 
which is in accordance with ()3.25p . 

Assemblies Let Zj ~ Po{aj), aj > 0, j > 1 . In this case, 



and 



«? = tT' S(^'^(a;) =exp(a,x^), j>l 



I 



g{x) = exp j , g'^^Kx) = exp j , ^ > 1, = aj, j > 1. 

Thus, the radius of convergence of equals to 00 for all finite / > 1. Assuming g G RTp, < p < cxd, 
we then apply Schur's lemma to (j4.34p with / = gf^i-, , fi = g and /2 = to get from the definition 
p.25p that < g^'-' < 00. By virtue of Theorem 1, this says that 



Proposition 2 An assembly converges if and only if the sequence {aj, j > 1} is such that g G RTp, 
for some < p < 00. 



Corollary 2 Assemblies with the following parameter functions a = {aj, j > 1} are convergent: 
(i) Smoothly growing: a £ RTp, < p < 00; 

(a) Oscillating: a G ^i^' + ^^"^'^v)^ where r>0, 0<e<| and y > 0. 

Proof (i) follows from the important for our study Corollary 4.3 in [6], which says that a £ RTp, < 
p < 00 implies g £ RTp, with the same p. For the proof of (ii), we employ the asymptotic formula 
(4.99) in [19] which says that for the assembly in question with y = 1 we have g G RTi, and then 
apply (ii) of Remark 2 with 9 = y. ■ 
Note that the restriction on e in the part (ii) of the last corollary determines a bound on the 
"size" of oscillation of TEZj = aj, j > 1 that ensures the RTp property for g. 

Examples. In combinatorics (see Table 2.2 in [3]), many assemblies, e.g. permutations (oj = i), 
Ewens sampling formula (aj = j, > 0), forests of labelled rooted trees {aj = ~ const e^j'^^), 
etc are regularly varying, which implies that that in these cases a G RTi, y > 0, where y is as in 

y 

the definition of a regularly varying structure. Note that in statistical mechanics, regularly varying 
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assemblies with y = 1 and a > are called generalized Maxwell-Boltzman statistics ([33|- Thus, by 
virtue of the condition (i) of Corollary 2, all the aforementioned assemblies converge. 

We give now two examples of divergent assemblies. Firstly, set partitions {aj = J — 1) 

diverges, since in this case g{x) = e^"^"^, so that the radius of convergence of g is infinity. In this 
connection note that the following sufficient condition for RT^q property of g was recently established 
in [TT]: If the parameter function a is such that gcd{j : aj > 0} = I and aj = 0{j^^ /jl), < 6 < 
1, j ^ oo, then g £ RToo, which means that the induced assembly diverges. 

For our second example we construct an assembly with g that does not belong to any class RTp, < 
p < oo. It is clear that the corresponding aj, j > 1 should exhibit a wild behavior. We set 

if i > 1 is odd, 
j-i + 2J+^j~\ if J > 2 is even. 



1 



1 - 4x2 



4I+K2] _ I 



We then have 

which by Cauchy product formula gives 

[n/2] 
k=0 

Consequently, 

lini = 1/4 

lim '^^^ = 1. 

C2„+l 
(2) 

Finally, note that for graphs on n vertices, c„ = ^-fj-, so that Cn € RTq and by (i) of Remark 1, for 
any / > 1 the limit measure is concentrated on the singleton (0, . . . , 0) G M;. 

We shift now to considering CFP's related to assemblies. In this case the ratio of the net 
transitions (jl.lip has the following form: 



where 



if i^j: hk^ > 

aiaj ' ^ J J J (4.37) 

^ if i = i : fei > 2, 



1: 



if i ^ j '■ kikj > 



V{ki,kj,ki+j) = { ^ ' ^ ' ' (4.38) 

^fe^, if i = 3^h>'2. 

As we explained in Section 1, such ratios correspond to mean field CFP's with net transition rates 

of coagulation and fragmentation as given by (I1.12|) . ()1.13p . Consequently, the preceding discussion 

reveals that amongst mean field CFP's both convergent and divergent models exist. 
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Multisets Assuming Zj ~ NB{'mj,p'), < p < 1, mj > 1, j > 1, we have 



where p is a free parameter. This gives 

4'^ = '^^ " s^'H^) = (1 - \xp\ < 1, 

which leads to the Euler type generating function 

~g{x) = ll{i-{pxyy"'\ (4.39) 

It is known (see e.g. |10J, Lemma 1.15) that g converges at some x : \xp\ < 1 if and only if 

lim (i-ipxyy""' = 1. 

This tells us that g has the same radius of convergence p as the series 

rrijpPxK (4.40) 

Since it is assumed that rrij > 1, i > 1, we have < p < p^^. Next, in ()4.34p the function 
= JiT(x) = nj=i (1 - {pxyy^' , / > 1 is a polynomial, such that /2(p~^) = 0, while 72(2;) > 
0, < X < p-^. 



Proposition 3 A multiset is convergent if and only if g £ RTp, for some p < p ^. 

Proof. In the case g G RTp-i, we have g^'^ = 0, Z > 1, by (|4.36|) and the observation following 
(j4.40p . So, in this case a multiset diverges, because condition (c) of Theorem 1 is violated. If now 
g G RTp with < p < p^^, then < qi < 00, / > 1, by the same argument. ■ 

Corollary 3 Multisets with the following parameter functions m = {mj, j > 1} are divergent: 

(i) m ^ RTi; (ii) mj = 0{j°'), for some a ^ R; (Hi) m G + ^5 !)> < e < |, r > 0, 

whereas multisets with parameter functions (iv) and (v) below converge: 

(iv) m G RTp, for some < p < 1; (v) m & -^(t + ^1 '''i 2/)' < e < |, y > 1, r > 0; 

(vi) rrij X y^j^^,y > 1. 

Proof. Each one of the conditions (i)-(v) is sufficient for g G RTp with a corresponding p. Namely, the 
conditions (i) and (iv) are due to Bell-Burris Theorems 6.1 and 6.2 ([6]), which state that m G RTp, 
for some < p < 1 implies g{p~^x) G RTp with the same p, and consequently, g{x) G RTpp-i. (Note 
that in Theorem 6.1 in [6], condition (c) is required only for the second part of the claim). 
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The condition (ii) provides g{p~^x) G RTi, by the powerful result of Bell [5] that generalizes 
Bateman and Erdos theorem. So, under this condition g(x) £ RTp-i. The conditions (iii) and (v) 
result from Corollary 2 of [23j for expansive multisets, which says that in the both cases g{p~^x) £ 
RTy-i, which is equivalent to g{x) S RTy-ip-i, with y > 1. Regarding the condition (vi), we firstly 
recall that aj x bj means that the ratio j > 1 is bounded above and below by positive constants. 
The sufficiency of the condition (vi) for g G RTy-ip-i, y > 1 was proven in [35] devoted to logical 
limit laws for logarithmic multisets. ■ 

Examples. Integer partitions {nij = 1, j > 1), planar partitions {nij = j, j > 1) (see [Ij) and 
generalized Bose- Einstein statistics {nij = j", a > 0) diverge, since in all these cases m G RTi. 
The following logarithmic multisets (see [3]) with m £ RTp, p < 1 converge : mapping patterns 
(nij = b < 1, j > 1), monic polynomials over GF(q), q > 1 {rrij ~ y). Also converge forests of 
unlabeled unrooted {mj ~ const b'^j"^^"^, b < 1) and rooted {mj ~ const b^^j^^^'^, b < 1) trees. 

Remark 3 There is a formal linkage between assemblies and multisets, expressed as follows. We 
saw that in the case of assemblies the generating fuction g is of the exponential form: g{x) = e^^^\ 
where Q{x) = X^j>i ajX^ , with aj > 0, j > 1. For multisets, the Euler type generating function g in 
can be written in the same form, with the function Q* called the star transformation (see f^, 
J^) of the generic generating function Q{x) = Sj>i mjp'x^ of the sequence m : 

Q*{x)=^m*x\ ^j = ^—f-^ J > 1, mg = 0. 

j>i ik=j 

Thus, m* > 0, J > 1, which is a basic assumption in the theory of RTp sequences. It was proven in 
that if mjp^ S RTp, with some < p < 1, then ~ rujp^ , j oo, which means that in this 
case the multiset behaves asymptotically as the assembly induced by Zj ~ Po{mjp'), j > 1. This fact 
explains the condition (iv) of Corollary 3. 

Regarding the associated CFP's, for multisets (jl.lip becomes: 

q{r^-t,j)=V{h,kj,h+,) I if ^/J :/ci/c, >0 ^^^^^^ 

where V{ki, kj, k-i^j) is as in (j4.37p . Since the second factor in (j4.4ip depends both on rj — (^1; • • • 5 ^n) 
and the parameters mj, j > 1 of the CFP, the process is not a mean field model. As we mentioned 
in section 1, this fact is especially clear in the case mj = 1, j > 1 that corresponds to the uniform 
measure /i„ on the set so that we have from (j4.4ip 

Qimhj) = 1, e ^n, i,j i + j< n. 
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Selections. In this case Zj ~ Bi{j^^;mj), rrij > 1, j > 1 and p > 0. Hence, 

So, as in the case of multisets, g has the same radius of convergence p < p~^ as the series (j4.40p . 
while the radius of convergence of g^^^ equals oo for all / > 1. 



Proposition 4 A selection is convergent if and only if the sequence {'rnj}j>i is such that g G RTp, 
for some < p < p^^. 

Proof: Since the radius of convergence of is p~^, for alH > 1, the preceding discussion implies 
that Schur's lemma can not be applied to (|4.34p . So, we apply the lemma to (|4.33p to get 

oo > g(') = (^('^(p)) ^>0, 1>1, 0<p<p-^. 



Now the claim follows from Theorem 1. ■ 
Comparing Propositions 3 and 4, one sees that for selections the conditions of convergence are less 
strict than the ones for multisets. 

For selections, conditions on mj, j > 1 of Bell-Burris type providing g G RTp are not known. 
This is due to the fact that for selections the star transformation does not ensure the nonnegativity 
of m*, j > 1. 

However, the probabilistic method for enumeration of decomposable structures works in the 
case of selections either. It follows from Theorem 5 in |23] that for expansive selections, m G 
+ e,r;y), r,y > 0, implies that g G RTp, with p = (yp)^^- By Proposition 4, this provides 
convergence of the corresponding selections, if y > 1. As a result, we derive the convergence of the 
following expansive selections: integer partitions into distinct parts {mj = 1, j > 1) and generalized 
Fermi statistics {mj = j", a > 0). 

For the associated CFP's we obtain from (jl.lip 



+1 



-pr, if i 7^ 7 : hkj > 0, < h < 



m2i—k2i 



{mi-ki+l){mi-ki+2) 

This shows that these CFP's are not mean field models. 



if i = j : 2 < ki < mi 



5 Concluding remarks and history 

Generally speaking, the phenomenon of asymptotic independence of a finite number of small groups 
of particles in large random systems (i.e. systems formed of a large number of randomly interacting 
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particles) was observed in different fields of applications, under various mathematical settings. The 
assumption of asymptotic independence, sometimes accepted without proof, was of great help for 
the study of probabilistic models considered. Without surveying the subject, we point out below a 
few settings parallel (in some sense) to the one in the present paper. 

(i) Mean field models. In statistical physics, there is a general belief that when the range of 
interaction between individual groups of particles is large enough, the aforementioned phenomenon 
of asymptotic independence occurs. The assumption of independence is the main idea behind the 
Ginzburg-Landau equations and it is the basis of calculations in a variety of mean field models and 
models that can be approximated by a mean field (see e.g. \T^, p. 111-113). Recall that in mean 
field models, each group of particles is allowed to interact with all other ones, so that the topology 
of a state space does not play a role. An example of a mean field is an interacting spin system on a 
complete graph with n vertices. Such a system is given by a time continuous Markov process on a 
state space {0, 1}", with given infinitesimal (in time) rates of flips — > 1 and 1 — > (for references 
on interacting particle systems see [28j). It is natural to suggest that under some conditions on 
the rates of flips, the equilibrium states of a given finite set of vertices of a complete graph become 
asymptotically independent, as n ^ oo. For a spin system called the noisy voter model the above 
fact was proven in [22] for pairs of states. (The pairwise independence of states at equilibrium follows 
from formula for covariance in Example 1, p. 38 in [22], after adding the missing factor r = — 1 in 
the exponent of the integrand and passing to the limit as N,t ^ oo). A good discussion of mean 
field approximations for systems of interacting particle is presented in [9]. 

From the above description, it is clear that in contrast to the setting (jl.ip in the present paper, the 
formulation of a mean field model does not assume any conditioning. 

While mean field theory arose in statistical physics it has been recently applied elsewhere, for 
example in neural networks, fluid mechanics and artificial intelligence, where it is used as an ap- 
proximation of complex (in respect of interactions) many-body systems in equilibrium. The main 
idea of the mean field approximation is to replace all interactions to any one body with one average 
interaction. As a result, the mean field partition function associated with a hamiltonian of a system 
considered becomes a product of the mean field local partition functions, which greatly simplifies the 
calculations (for references see e.g. Section 6.1.4 in |29j). 

(ii)The Gibbs conditioning principle (see [16], [13], [E]). In the context of an ideal gas model, 
the simplest version of the principle reads as follows. Let Xi,X2, ... be independent and identically 
distributed random variables viewed as energies of individual particles, so that Xi + . . . X^ is the 
total energy of a system of n particles. Let lEXi = 1 and assume some suitable regularity conditions 
on a common probability law P of the sequence of random variables. Then, for a fixed k > \ the 
distribution law Afc^„ of (Xi, . . . , Xk) conditioned on Xi -|- . . . Xn = n weakly converges, as n ^ oo to 
the k— fold product law . In statistical physics the law A^ „ is called a microcanonical distribution. 
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and the gibbs conditional principle asserts asymptotic independence of energies of any finite number 
of particles in microcanonical ensembles. Formulated in the beginning of the 20-th century, the 
principle has been extended and refined in different directions, with particular attention being paid 
to the rate of convergence to limit distributions. 

To distinguish from the conditioning relation (jl.ip . the measure ^k,n is defined on the simplex, 
from which it follows that the Cauchy product relationship (|4.33p . which is basic for the study of 
multiplicative measures is not valid. However, we believe that the interplay between the above 
two settings deserves further study. Quite independently, the distribution Afc „, with k = n and 
discrete and not necessary identically distributed random variables Xi , X2 , ■ ■ ■ was introduced by 
Kolchin ([27]) for representation of distributions of cell counts in combinatorial urn schemes. In 
|27j the representation is called the generalized scheme of allocation, whereas in [7J and [30j, it is 
named Kolchin representation formula. We note that the problem of asymptotic independence is not 
addressed in [27] . 

(iii) Random combinatorial structures. In the theory of random structures, the asymptotic 
independence of numbers of small groups was discussed in numerous papers, starting from the 1940's. 
A general set up leading to asymptotic estimation of the total variation distance between component 
spectrum of small counts (as defined by the conditioning relation (jl.ip ) and the independent process 
was developed by Arratia and Tavare in their seminal paper [2] (see also [3], [21] and [27]). As a 
result, asymptotic independence was established for logarithmic random structures with y = 1. In 
[1], the same was proven for convergent structures, and in [19], the asymptotic independence was 
proven for expansive assemblies with any y > 0. Regarding assemblies, multisets and expansive 
selections, the aforementioned results easily follow from our results in Section 4. In fact, recall that 
for assemblies EZj = aj, j > 1„ for multisets lEZj = jz^, i^l) 0<p<l and for selections 
TEZj = j >1, < p < 1. By the definition (OHD of a regularly varying structure, we thereby 

conclude that the following facts hold. 

• For assemblies, a = {aj, j > 1} £ RTy-i, y > 0, for all a G IR, which implies convergence by 
the condition (i) of Corollary 2. 

• For multisets, m = {rrij, j > 1} G RTy~ip, < p < 1, for all a £ R. This implies convergence 
if y ^ 1) by the condition (iv) of Corollary 3. Note that in the case y = p, a regularly varying 
multiset diverges for all a G H, by the condition (i) of Corollary 3. 

• For selections, nij ~ constj°'{yp^^y , j > 1. By the discussion following Proposition 4, this 
provides convergence for all a > and y > 1. 

In this connection, we mention that to our knowledge no examples of random structures for which 
the independence principles fails were given in the literature, prior to this paper. 
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Our results reveal also a basic difference between pictures of asymptotics clustering of components 
in general convergent structures and regularly varying divergent multisets. Namely, in the case of a 
convergent structure, Corollary 1 tells us that with a positive limit probability there are components 
of any fixed sizes, i.e. 

lim Ffiicj"^ =A;i,...,k/"^ = A:/) >0, ki,...,ki>0, I > 1. 

n—>oo 

On the other hand, if a regularly varying multiset diverges, then g*-'^ =0, / > 1, by Proposition 3, 
from which it follows that lim^^oo ]P(-^j"'* = kj) = 0, for all finite kj > 0, j > 1. 

(iv) CFP's. It is common to trace the beginning of rigorous mathematical models of coagulation- 
fragmentation to the paper by Smoluchowski (1918) where the famous system of coagulation equa- 
tions describing the time evolution of the process was derived. Already in this paper the assumption 
of independence (more precisely, the absence of correlations) of clusters of small sizes was adopted. 
Subsequently, deterministic and stochastic versions of the model were studied in numerous papers 
in probability and various applied fields. The study of reversible CFP's was concentrated on what 
we call in the present paper mean field CFP's (|35j. |26j). (Recall that these models conform to 
assemblies) . In jl5j the model was treated as a reversible Markov chain on the set of partitions and 
it was proven (Theorem 4, (4.24)) that, if ^ G RTp, < p < oo, then at the equilibrium of the 
process, 

cov{kI''\ K^^) ^0, n ^ oo, 

for any fixed / 7^ m. This is, of course, a weak form of our Proposition 2. More details on the history 
of CFP's can be found in [19j. 
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